H,O, Sr,Nd, and Pb isotope geochemistry of the Latir volcanic field and cogenetic intrusions,New Mexico,and relations between evolution of a continental magmatic center and modifications of the lithosphere

Over 200 H, O, Sr, Nd, and Pb isotope analyses, in addition to geologic and petrologic constraints, document the magmatic evolution of the 28.5–19 Ma Latir volcanic field and associated intrusive rocks, which includes multiple stages of crustal assimilation, magma mixing, protracted crystallization, and open- and closed-system evolution in the upper crust. In contrast to data from younger volcanic centers in northern New Mexico, relatively low and restricted primary ¹⁸O values (+6.4 to +7.4) rule out assimilation of supracrustal rocks enriched in ¹⁸O. Initial ⁸⁷Sr/⁸⁶Sr ratios (0.705 to 0.708), ¹⁸O values (-2 to-7), and ²⁰⁶Pb/²⁰⁴Pb ratios (17.5 to 18.4) of metaluminous precaldera volcanic rocks and postcaldera plutonic rocks suggest that most Latir rocks were generated by fractional crystallization of substantial volumes of mantle-derived basaltic magma that had near-chondritic Nd isotope ratios, accompanied by assimilation of crustal material in two main stages: 1) assimilation of non-radiogenic lower crust, followed by 2) assimilation of middle and upper crust by inter-mediate-composition magmas that had been contaminated during the first stage. Magmatic evolution in the upper crust peaked with eruption of the peralkaline Amalia Tuff (26 Ma), which evolved from metaluminous parental magmas. A third stage of late, roofward assimilation of Proterozoic rocks in the Amalia Tuff magma is indicated by trends in initial ⁸⁷Sr/⁸⁶Sr and ²⁰⁶Pb/²⁰⁴Pb ratios from 0.7057 to 0.7098 and 19.5 to 18.8, respectively, toward the top of the pre-eruptive magma chamber. Highly evolved postcaldera plutons are generally fine grained and are zoned in initial ⁸⁷Sr/⁸⁶Sr and ²⁰⁶Pb/²⁰⁴Pb ratios, varying from 0.705 to 0.709 and 17.8 to 18.6, respectively. In contrast, the coarser-grained Cabresto Lake (25 Ma) and Rio Hondo (21 Ma) plutons have relatively homogeneous initial ⁸⁷Sr/⁸⁶Sr and ²⁰⁶Pb/²⁰⁴Pb ratios of approximately 0.7053 and 17.94 and 17.55, respectively. ¹⁸O values for all the postcaldera plutons overlap those of the precaldera rocks and Amalia Tuff, except for those for two late-stage rhyolite dikes associated with the Rio Hondo pluton that have ¹⁸O values of-8.6 and-9.5; these dikes are the only Latir rocks which may be largely crustal melts.Chemical and isotopic data from the Latir field suggest that large fluxes of mantle-derived basaltic magma are necessary for developing and sustaining large-volume volcanic centers. Development of a detailed model suggests that 6–15 km of new crust may have been added beneath the volcanic center; such an addition may result in significant changes in the chemical and Sr and Nd isotopic compositions of the crust, although Pb isotope ratios will remain relatively unchanged. If accompanied by assimilation, crystallization of pooled basaltic magma near the MOHO may produce substantial cumulates beneath the MOHO that generate large changes in the isotopic composition of the upper mantle. The Latir field may be similar to other large-volume, long-lived intracratonal volcanic fields that fundamentally owe their origins to extensive injection of basaltic magma into the lower parts of their magmatic systems. Such fields may overlie areas of significant crustal growth and hybridization.     

Implications of abundant hygroscopic minerals in the Martian regolith

Converging lines of evidence suggest that a significant portion of the Martian surface fines may consist of salts and smectite clays. Salts can form stoichiometric hydrates as well as eutectic solutions with depressed freezing points; clays contain bound water of constitution and adsorb significant quantities of water from the vapor phase. The formation of ice may be suppressed by these minerals in some regions on Mars, and their presence in abundance would imply important consequences for atmospheric and geologic processes and the prospects for exobiology.     

Martian glaciation and the flow of solid CO2

The flow law determined experimentally for solid CO₂ establishes that a hypothesis of glacial flow of CO₂ at the Martian poles is not physically unrealistic. Compression experiments carried out under 1 atm pressure and constant strain rate demonstrate that the strength of CO₂ near its sublimation point is considerably less than the strength of water ice near its melting point. The data fit a power law “creep” equation of the form

$\text{?}\text{?}\text{= (4 × 10}{}^6\text{) σ}{}^{3.9}\text{exp}\text{(}\text{?12 200}\text{RT}\text{)}$

, where ? is compressive strain rate (sec^?1), σ is compressive stress (bars), R is the gas constant in calories per mole, and T is absolute temperature. The exponent of σ of 3.9 contrasts with a value near 3.1 for water ice, and indicates that the strain rate is somewhat more sensitive to stress for CO₂ than for water. Likewise, the low activation energy for creep, 12 200 cal mole^?1, illustrates that CO₂ is not highly sensitive to temperature and is thus likely to flow over a broad range of temperatures below its melting point. Strength values for CO₂ are of the order of one-tenth to one-third the strength of ice under equivalent conditions.A plausible glacial model for the Martian polar caps can be constructed and is helpful in explaining the unique character of the polar regions. CO₂-rich layers deposited near the pole would have flowed outward laterally to relieve high internal shear stresses. The topography of the polar caps, the uniform layering of the layered deposits, and the general extent of the polar “sediments” could all be explained using this model. Flow of CO₂ rather than water ice greatly reduces the problems with Martian glaciation. Nevertheless, problems do remain, in particular the large amounts of CO₂ necessary, the need to increase vapor pressure and temperature with depth in the polar deposits, and the lack of good observational evidence of flor features. Within the limits of the present knowledge of surface conditions of Mars, CO₂ glaciation appears to be a realistic alternate working hypothesis for the origin of the polar features.     

Surface properties of asteroids: A synthesis of polarimetry,radiometry, and spectrophotometry

The surface compositions of 110 asteroids are analyzed from statistically representative data sets of polarimetry as a function of phase angle, broad-band radiometry near 10 and 20 μm, and visible and near-infrared spectrophotometry. A comparison of albedos and diameters determined by polarimetry and radiometry shows that a modest upward revision of the radiometric albedo scale is needed and that a single law relating the slope of the polarization-phase curve to geometric albedo may not hold for very dark asteroids. We present reliable adopted albedos and diameters for 56 objects. Roughdi ameters for 52 additional objects are obtained from spectrophotometry using a correlation between albedo and color. Corrections for sampling bias permit investigation of asteroid compositions as a function of diameter, orbit, and other parameters.More than 90% of the minor planets fall into two broad compositional groups, defined by several optical parameters, designated by the symbols C and S. Comparisons with meteorite spectral albedo curves suggest that the two groups are compositionally similar to carbonaceous and stony-metallic meteorites, respectively. C-type asteroids predominate in the belt, especially in the outer half. An unusual distribution of compositions is found between 2.77 and 3.0 AU. Many S-type objects have diameters of 100–200 km; C-type objects are much more common at both larger and smaller sizes. Vesta is unique, being apparently the only differentiated asteroid remaining intact in the belt. The largest C-type objects are compositionally distinct from smaller ones and possibly are metamorphosed. We sketch some implications for meteoritics and for the early history of the solar system and point to the need for further systematic sampling of smaller and fainter objects by these three observational techniques.     

The spectral reflectivity of select areas on Europa

The July 12, 1973, occultation of Europa by Io was observed in 30 wavelength channels in the spectral region λλ3200-11, 000 Å with the 200-inch Hale telescope and a multichannel spectrometer. The data are presented in absolute units above the Earth's atmosphere. The data are analyzed to obtain the spectral reflectivity of seven localized areas on the disk of Europa centered on 324°W longitude. The equatorial material is confirmed to be darker than the eastward-skewed bright north polar cap and a hint is evident that the darker material as well may be somewhat redder than the cap material.     

Diffusion-limited REE uptake by eclogite garnets and its consequences for Lu–Hf and Sm–Nd geochronology

Garnets from the Zermatt-Saas Fee eclogites contain narrow central peaks for Lu + Yb + Tm ± Er and at least one additional small peak towards the rim. The REE Sm + Eu + Gd + Tb ± Dy are depleted in the cores but show one prominent peak close to the rim. These patterns cannot be modeled using Rayleigh fractionation accompanied by mineral breakdown reactions. Instead, the patterns are well explained using a transient matrix diffusion model where REE uptake is limited by diffusion in the matrix surrounding the porphyroblast. Observed profiles are well matched if a roughly linear radius growth rate is used. The secondary peaks in the garnet profiles are interpreted to reflect thermally activated diffusion due to temperature increase during prograde metamorphism. The model predicts anomalously low ¹⁷⁶Lu/¹⁷⁷Hf and ¹⁴⁷Sm/¹⁴⁴Nd ratios in garnets where growth rates are fast compared to diffusion of the REE, and these results have important implications for Lu–Hf and Sm–Nd geochronology using garnet.     

Seismic Hazards Assessment for Radioactive Waste Disposal Sites in Regions of Low Seismic Activity

A comparative tectonic quiescence and lack of earthquakes make the stable centres of continents attractive for siting long-term radioactive waste storage facilities. The low rates of deformation in such regions, however, make it difficult to characterize their long-term seismotectonic behaviour, leading to uncertain estimates for the very low probability hazard estimates required by society. In an attempt to overcome the deficiency of both contemporary seismicity and paleoseismic data in central Canada, we have used earthquake histories from regions with similar seismotectonic characteristics from around the world. Substituting space for time, we estimate a long-term rate per 10⁶ km² of 0.004 magnitude ≥6 earthquakes per annum, of which 33–100% might rupture to the surface.     

Groundwater residence times and flow paths in fractured rock determined using environmental tracers in the Mission Tunnel; Santa Barbara County,California, USA

Multiple geochemical tracers [ion chemistry, stable isotopes of water, chlorofluorocarbons (CFC), tritium] and a 25-year-long record of discharge were used to understand residence times and flow paths of groundwater seeps in the fractured rock aquifer surrounding the Mission Tunnel, Santa Barbara, California. Tritium data from individual seeps indicate that seep waters are a mixture of >45-year-old (recharged prior to the nuclear bomb tests) and young groundwater. CFC data support this interpretation, however, a two-end member mixing model cannot completely explain the age tracer data. Microbial degradation and partial re-equilibration complicate the CFC signal. Spectral analysis of precipitation and groundwater seepage records shows that seepage lags precipitation by 3 months. This delay is related to the advancement of the wetting front and increasing the number of active flow paths. Additionally, the amount of seepage produced by precipitation is less during extended periods of drought than during normal or wet periods, suggesting antecedent conditions strongly affect flow through this fractured rock aquifer.     

Methods for the visualization of digital elevation models for landform mapping

Digital elevation models (DEMs) are increasingly used for landform mapping, particularly with the growing availability of national and global datasets. In this paper we describe a variety of techniques that can visualize a DEM. We then compare five techniques to ascertain which performs the most complete and unbiased visualization. We assess the visualization techniques by comparing landforms mapped from them against a detailed morphological map (derived from mapping of multi‐azimuth relief‐shaded DEMs cross‐checked with stereo aerial photographs). Results show that no single visualization method provides complete and unbiased mapping. The relief‐shaded visualizations are particularly prone to azimuth biasing, although they can highlight subtle landforms. We recommend curvature visualization for initial mapping as this provides a non‐illuminated (and therefore unbiased) image. Initial mapping can then be supplemented with data from relief‐shaded visualizations. Copyright © 2005 John Wiley & Sons, Ltd.     

Multi‐year water balance assessment of a newly constructed wetland,Fort McMurray,Alberta

Oil sands mining in Alberta transforms the boreal landscape of forests and wetlands into open pits, tailings ponds and overburden piles. Whereas reclamation efforts have primarily focused on upland forests, rebuilding wetland systems has recently become a motivation for research. Wetland creation and sustainability in this region is complicated by the sub‐humid climate and salinity of underlying mining material. In 2012, Syncrude Canada Ltd. completed the construction of the Sandhill Fen Watershed (SFW), a 52‐ha upland‐wetland system to evaluate wetland reclamation strategies on soft tailings. SFW includes an active pumping system, upland hummocks, a fen wetland and underdrains. To evaluate the influence of management practices on the hydrology of the system, this study reports the water balance from January 2013 to December 2014, the first 2 years after commissioning. A semi‐distributed approach was taken to examine the fluxes and stores of water in uplands and lowlands. Natural and artificial inputs and outputs were measured using a series of precipitation gauges and pumps, and evapotranspiration was quantified using three eddy covariance towers. A series of near surface wells recorded water table position. Both 2013 and 2014 were normal rainfall years, with 2013 having more and 2014 less snow than normal. In 2013, inflow/outflow from pumping was the predominant hydrological fluxes, resulting in considerable variability in water table position and storage changes throughout the summer. In 2014, the artificial addition of water was negligible, yet the water table remained near the surface in lowland locations, suggesting that wetland conditions could be maintained under current conditions. Evapotranspiration rates between uplands and lowlands were similar between years and sites, ranging from 2.2 ± 1.8 to 2.5 ± 1.2 mm/day and were largely controlled by climate. These rates were less than nearby older upland systems, suggesting that water balance partitioning will change as vegetation develops. Comparison between years and with natural systems provides insight on how management practices influence hydrologic dynamics and the overall water balance of the SFW. Copyright © 2016 John Wiley & Sons, Ltd.